Has the QCD Critical Point been Signaled by Observations at RHIC ? 
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The shear viscosity to entropy ratio (rj/s) is estimated for the hot and dense QCD matter created 
in Au+Au collisions at RHIC {^/snn = 200 GeV). A very low value is found rj/s ~ 0.1, which 
is close to the conjectured lower bound (l/47r). It is argued that such a low value is indicative of 
thermodynamic trajectories for the decaying matter which lie close to the QCD critical end point. 
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A primary objective for studying ultra-relativistic 
heavy ion coUisions is to map out the accessible regions 
of the Quantum Chromo Dynamics (QCD) phase dia- 
gram. The existence of a critical point - an end point of 
the first order chiral transition in this diagram - has at- 
tracted considerable attention since it was first proposed 
[1, 2]. Several theoretical techniques have been exploited 
to locate this critical end point (CEP) in the plane of 
temperature vs baryon chemical potential (T, /is) i.e. the 
QCD phase diagram [3, 4, 5, 6, 7] . The properties of the 
CEP have also been sought via model calculations and 
the universality hypothesis [8, 9, 10]. Nonetheless, ex- 
perimental verification of the CEP constitutes a major 
current scientific challenge. 

The study of heavy ion collisions to search for the 
CEP was proposed several years ago [11]. Recently, a 
resurgence of experimental interest has centered around 
the prospects for successful experimental searches at 
Brookhaven's Relativistic Heavy Ion Collider (RHIC) 
[12, 13] and CERN's Super Proton Synchrotron (SPS) 
[14, 15]. Indeed, a recent theoretical investigation in- 
volving two flavors of light dynamical quarks [IG] gives 
Tc/Tco 0.95 and n%/Tco ~ 1.1, where ~ 170 McV 
is the temperature for a cross-over to the quark gluon 
plasma (QGP) at zero baryon chemical potential and Tc 
and fig are the temperature and the baryon chemical 
potential (respectively) at the critical end point. It is 
widely believed that these values for Tc and fig place the 
CEP in the range of observability for future energy scans 
at both RHIC and the SPS [13]. 

In recent hydrodynamical calculations [17, 18], an in- 
teresting aspect of the CEP has been reported which 
could have important consequences for its detection. 
That is, it acts as an attractor for thermodynamic trajec- 
tories of evolving hot and dense QCD matter [11]. This 
aspect of the CEP is illustrated in Fig. 1 where isentropic 
trajectories in the (T, /xb) plane are shown for a broad 
range of values for the entropy per baryon s/hb- Trajec- 
tories for the range s/ub ~ 100 — 50, are clearly focused 
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FIG. 1: (Color online) Results from a hydrodynamical calcu- 
lation by Nonaka et al which assumes an EOS with a CEP 
[17]. Isentropic trajectories are shown for the case in which 
the CEP is located at (T,/is) = (154.7 MeV, 367.8 MeV). 
The value of s/ns for each of these trajectories (left to right) 
is 100.0, 66.6, 50.0, 40, 33.3, 28.6 and 25. 



toward the CEP. Such a focusing is not observed for isen- 
tropic trajectories in similar model calculations which do 
not take explicit account of the CEP in the equation of 
state (EOS). Instead, trajectories are merely shifted on 
the phase coexistence line [17]. Given this aspect of the 
CEP, it is worthwhile to investigate whether or not the 
decay dynamics of the high energy density QCD matter 
created in heavy ion RHIC collisions show any signals for 
the CEP. 

Here, we argue that estimates for the ratio of viscosity 
to entropy density (77/s) made from RHIC measurements 
[19, 20, 21], indicate decay trajectories which are close to 
the CEP. This implies that the prospects for constraining 
the location of the CEP in the (T, /is) plane are very 
good. Possibly, such constraints can be best achieved via 
a study of excitation functions with relatively large steps 
in collision energy prior to focusing on a series of fine 
steps. 
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FIG. 2: (Color online) vijeuq vs KEr/uq for several identi- 
fied particle species obtained in minimum bias Au+Au colli- 
sions. The data is taken from Refs. [32, 33] and the references 
therein. The dashed-dotted line represents a fit to the data. 



There is strong experimental evidence for the creation 
of locally equilibrated nuclear matter at unprecedented 
energy densities, in heavy ion collisions at RHIC [22, 23, 
24, 25, 26, 27, 28, 29]. Jet suppression studies indicate 
that the constituents of this matter interact with unex- 
pected strength and this matter is almost opaque to high 
energy partons [30, 31]. 

Elliptic flow measurements [32, 33] validate the pre- 
dictions of perfect fluid hydrodynamics for the scaling 
of the elliptic flow coefficient V2 with eccentricity e, sys- 
tem size and transverse kinetic energy KEt [34, 35, 3(3] ; 
they also indicate the predictions of valence quark num- 
ber ijiq) scaling [37, 38, 39], suggesting that quark-like 
degrees of freedom are pertinent when elliptic flow devel- 
ops [32, 33]. The result of such scaling is illustrated in 
Fig. 2 where we plot V2/nq£ vs KEt/jIij; it shows that the 
relatively "complicated" dependence of V2 on centrality, 
transverse momentum, particle type and quark number 
can be scaled to a single function [32, 33]. 

The expected suppression of flow due to the shear vis- 
cosity, as predicted by weak-coupling transport calcula- 
tions is also not observed. In fact, only a rather small 
viscosity, ??/s < 0.1, could be accommodated by the data 
in early estimates [19]. This value of 77/ s, which is much 
lower than that obtained from weak coupling QCD [40] or 
hadronic computations [41], has been interpreted as evi- 
dence that the quark gluon plasma (QGP) created in the 
early phase of RHIC collisions is more strongly coupled 
than expected [28, 42]. An alternative interpretation is 



that the low value for 77/s could result from an anoma- 
lous viscosity r|A^ arising from turbulent color magnetic 
and electric fields dynamically generated in the expand- 
ing quark-gluon plasma [43]. That is, ■q~^ = r]^^ Vc^ 
and rjA dominates over the coUisional viscosity rjc ■ 

The ratio rj/s cannot be arbitrarily small because 
quantum mechanics limits the size of cross sections via 
unitarity. A conjectured lower bound for r]/s is l/47r, 
reached in the strong coupling limit of certain gauge the- 
ories [44]. It has even been speculated that this lower 
bound holds for all substances [44] . The temperature de- 
pendence of rj/s also provides invaluable insights. This 
is illustrated in Fig. 3 where we have plotted ry/s vs 
(T — Tc)/Tc for molecular, atomic and nuclear matter. 
The data for He, N2 and H2O were obtained for their re- 
spective critical pressure. They are taken from Ref. [45] 
and the references therein. The calculated results shown 
for the meson-gas (for T < Tc) are obtained from chiral 
perturbation theory with free cross sections [46]. Those 
for the QGP (ie. for T > Tc) are from lattice QCD sim- 
ulations [47]. The value Tc ~ 170 MeV is taken from 
lattice QCD calculations [48]. 

Figure 3 illustrates the observation that for atomic and 
molecular substances, the ratio 77/s exhibits a minimum 
of comparable depth for isobars passing in the vicinity of 
the liquid-gas critical point [44, 45, 46]. When an isobar 
passes through the critical point (as shown in Fig. 3), the 
minimum forms a cusp at Td when it passes below the 
critical point, the minimum is found at a temperature 
below Tc (liquid side) but is accompanied by a discontin- 
uous change across the phase transition. For an isobar 
passing above the critical point, a less pronounced min- 
imum is found at a value slightly above Tc- The value 
77/s is smallest in the vicinity of Tc because this corre- 
sponds to the most difficult condition for the transport 
of momentum [45] . 

Given the observations for atomic and molecular sub- 
stances, one expects a broad range of trajectories in the 
(T, ns) or (T, /is) plane for nuclear matter, to show rj/s 
minima with a possible cusp at the critical point. The 
exact location of this point ( in the (T, /xb) plane) is of 
course not known, and only coarse estimates of where 
it might lie are available. The open triangles in Fig. 3 
show calculated values for 77/s along the hb = 0, ns ~ 
trajectory. For T < Tc the 77/s values for the meson- 
gas show an increase for decreasing values of T. For T 
greater than Tc, the lattice results [47] indicate an in- 
crease of 7;/s with T, albeit with large error bars. These 
trends suggest a minimum for 77/s. 

The open hexagons in Fig. 3 show estimates for 7;/s 
constrained by experimental data; they are for a tra- 
jectory which starts along ub = nQ. The 77 values are 
extracted from momentum transport data obtained in 
intermediate-energy (IE) collisions [49, 50, 51]; values 
for s were obtained from composite yields measured in 
the same energy and mass region. The filled hexagon in 
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FIG. 3: (Color online) -q/s vs (T - T^)/Tc for several sub- 
stances as indicated. The calculated values for the meson-gas 
have an associated error of ~ 50% [4(i]. The lattice QCD 
value Tc = 170 MeV [ In] is assumed for nuclear matter. The 
lines are drawn to guide the eye. 

Fig. 3 shows a recent estimate for rj/s from high pre- 
cision elliptic flow measurements in Au-|-Au collisions 
(Vs^ = 200 GeV) at RHIC [21]: 

- ^ 1 XjCs, 

s 

where T is the temperature, A/ is the mean free path and 
Cs is the sound speed in the matter. 

The temperature T = 160 MeV was employed in a 
blast wave model generalized to include viscous correc- 
tions [19]. We estimate the value T = 165 ± 3 MeV via a 
hydrodynamically inspired fit to the universally scaled V2 
data shown in Fig. 2. The functional form used for the fit 
is /i {w) /lo (w) , where w = KEt /2T and Ii {w) and /q {w) 
are Bessel functions [35] . For Cs , we use the recent obser- 
vation that V2 scales with eccentricity and across system 
size to constrain the estimate Cs = 0.35±0.05 [32, 33, 52]. 
The mean free path estimate A/ = 0.3 ± 0.03/m [53], is 
obtained from an on-shell transport model simulation of 
the gluon evolution (in space and time) in Au-|-Au col- 
lisions at 200 GeV [54]. This parton cascade model in- 
cludes pQCD 2 2 and 2 <-> 3 scatterings. These values 
for A/, T and Cs lead to the estimate 77/s ~ 0.09 ± 0.015 
with an estimated systematic error of +0.1 primarily due 
to the uncertainty in A/. This value for rj/s is in good 
agreement with the experimentally based estimates of 



Teaney and Gavin [19, 20] and the theoretical estimates 
of Gyulassy and Shuryak [42, 55]. However, it is in stark 
contrast to the predictions of pertubative QCD [19, 40]. 

The experimental and calculated rj/s values for nuclear 
matter shown in Fig. 3, give strong indication for a min- 
imum in the vicinity of Tc- This minimum, achieved in 
RHIC collisions, is rather close to the conjectured lower 
bound oi rj/s = 1/4tt. We therefore conclude that it is 
compatible with the minimum expected if the hot and 
dense QCD matter produced in RHIC collisions follow 
decay trajectories which are close to the CEP. As dis- 
cussed earlier, such trajectories could be readily followed 
if the CEP acts as an attractor for thermodynamic tra- 
jectories for the decaying matter [17, 18]. Here, it is 
noteworthy that the system has to spend a considerable 
portion of its dynamics in the region of low j]/s for that 
low value to have a significant impact on the reaction 
dynamics. If the system gets momentarily, to a low rj/s 
and the rest of the dynamics takes place at high 77/ s, 
the momentary low value will have little impact. An op- 
timal situation would be if the system stays at low rj/s 
and then very quickly freezes out thereafter. That is, sig- 
nificant collective motion develops prior to the systems 
evolution toward higher rj/s values. 

In summary we have used the comprehensive overall 
scaling of elliptic flow data to constrain a rudimentary 
experimental estimate for 77/s in Au-|-Au collisions at 
RHIC [5()]. A value close to the conjectured lower bound 
oi rj/s = l/47r is observed. We argue that such a low 
value for 77/s is suggestive of decay trajectories which lie 
close to the QCD critical end point, because it acts as an 
attractor for constant entropy phase space trajectories 
of evolving hot and dense QCD matter. This suggests 
that the prospects for a successful search for the first 
order QCD phase transition and the possibility for es- 
tablishing reliable constraints (in the (T^^b) plane) for 
its associated critical end point are quite good. An initial 
measurement of excitation functions with coarse-grained 
collision energy steps may give sufficient information, or 
at least, give a better perspective on where to focus at- 
tention on fine steps. 
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